International Journal of Advanced Engineering Research and 
Science (IJAERS) 
ISSN: 2349-6495(P) | 2456-1908(0) 


msm hon 9$ Cras 
faia uw ND LM 


Y Vol-8, Issue-4; Apr, 2021 
riam aligai Jour mmu Journal Home Page Available: https://ijaers.com/ 
Engin iri ng Plaseuc* and Sheree T 
Journal DOI: 10.22161/ijaers 
Article DOI: https://dx.doi.org/10.22161/ijaers.84.17 





A study about Fourier series: Mathematical and graphical 
models and application in electric current and square 


oscillations 
Felipe Corrêa Veloso dos Santos!, Edivaldo Fernandes da Silva”, Anne Louise de Melo 
Dores”, Allan Thiago de Oliveira” 


-5^Department of Engineering, Pontifical Catholic University of Goiás, Brazil 
“Department of Mathematics UNU - Posse, University of the State of Goiás, Brazil. 





Received: 03 Jan 2021; Abstract— The present work aims to study the Fourier series with the 


pacsivad Inevised form intention of creating a didactic and understandable text. There are several 


ways to introduce the Fourier series, the simplest and most pedagogical 


12 Peb 2021; way, at least for beginners in the subject, is to explore intuitively the basic 
Accepted: 16 Mar 2021; idea through graphics. This analysis facilitates the understanding of the 
Available online: 16 Apr 2021 methodology behind the series and the Fourier. The mathematical 


©2021 The Author(s). Published by AI formulation will be a simple extension of the physical (geometric) image 


mE au l presented by the graphics. The popular saying goes that a figure is worth 
Publication. This is an open access article 


a thousand words. Mathematicians, in turn, argue that an equation is 
under the CC BY license worth a thousand figures. Maybe that is why in Fourier's analysis books, 
(https://creativecommons.org/licenses/by/4.0/). | mathematical formulation comes first, followed by the graphics. 

Keywords— Applications of Fourier 


series, Electric current, Fourier Series, 





Square oscillations. 





I. INTRODUCTION problem of heat conduction, which he explains in his work 
The theory of Fourier series began marked by the Théorie Analytique de la Chaleur in 1822, as the theory of 


attempt to solve heat conduction problems in a bar. In the heat diffusion in solid bodies. 


mid-eighteenth century a big problem arose that helped According to KOHAUPT and LEE (2015), the Fourier 
both mathematical analysis and other areas, including series was originated in telecommunication engineering, 
physics. The problem was to write a certain function given specifically in the connection with the digitalization of 
as an infinite series of sine and cosine functions, more data (signals) in order to transmit them. 


specifically a periodic function. Different from the current days where signals are 


Many physical phenomena related to celestial digitalized, formerly, signals were transmitted in the 
mechanics, particle mechanics, wave, ballistics, heat superposition of continuous periodic time functions, in the 
transmission, and electrical circuits are translated into form of waves. 


differential equations whose differential and integral Although Fourier’s proposition that any function 


calculus is not enough to solve them, arising the need for (continue or discontinuous) could be written as an infinite 


something more sophisticated, called Fourier equations. sum of the trigonometric functions cosine and sine had 


The development and creation of this series is due to been written down earlier, his contribution was ground- 
Jean Baptiste Joseph Fourier (1766-1830), French breaking by the degree of attention given to the 
physicist and mathematician. Fourier used it to solve the convergence of the functions (FOLKERTS et. al, 2020). 
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Still, according to AMORIM et. al (2019), although the 
work of Fourier had been heavily criticised by Lagrange, 
Laplace and Legendre, the fact that Fourier have observed 
that discontinuous functions are the sum of infinite series 1t 
was a great advance, for it explicit the shape of the series 
that represents this functions. 


The applications of Fourier series are the most varied. 
In the matter of electric current and oscillations, Fourier 
series are vital to approximate a periodic waveform in 
electronics and electrical circuits. For the modern 
technology and engineering, is by the use of this series that 
it is possible to decompose periodic signals into sum of 
infinite trigonometrically series in sine and cosine terms 
(ANUMAKA, 2012). 


In this work, we approach the Fourier series in three 
sections. In the first section, we deal with the subjects 
related to the definition of series, even and odd functions, 
periodic functions, as well as their properties. In the 
second, Fourier's trigonometric coefficients are 
determined, giving some examples of functions written in 
the form of the Fourier series. Finally, in the third section, 
we talk about some applications using odd and even 
function, and cases where the functions are neither even 


nor odd. 


II. METHODOLOGY 


In this chapter, we will study the definition of limits, 
derivatives, integrals, periodic functions, even functions, 
odd functions, sequences, series, their properties and their 
involvement with the series under study. We will now 
move on to formal settings. 


2.1 Limits 


We say that a function f (x) has limit ‘k’ when x — a, 
if for all e > O there is 6 > O as demonstrated by MUNEM 
and FOLIS (1882): 


0«lx-al«8 >| f co-kl «e (1) 


Such ‘k’ number when it exists is unique and will be 
indicated by: 


lim f (x) = k (2) 


2.1.1 | Limits properties 


For f (x) and g (x) functions of a real variable and K' 
belonging to IR. We consider: 


lim fix) =k (3) 
and 
lim g(x) =k (4) 
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Therefore: 


lim [f (x) + gtx] = lim flo + limg(x) = ky + ko: 
ii. lim [f (x) «glx)] = imf (x). limg(x) = ky * Ko: 

im [£2] — eh — desde que k + 0: 

eg "1 Siima k’ desde que Ka # O; 


i. lim[& = f()] 2 K&.limfix)- & ek. 





iii. 


All demonstrations here omitted can be seen in 
GUIDORIZZI (2001). 


2.2 Derivatives 


Among limits, there is a special called derivatives, 
which is our next object of study. Be f a function and k a 
point of your domain. The limit: 

lim fix-R)- f (x) 


EX 


(5) 


When it exists and is finite, it is named derivative and 
is indicated by f"(k). Therefore: 


nn — p Free fon 
f (X) = lim 


(6) 


If the limit above exists, we say that such function is 
differentiable at the given point. 


2.2.1] Derivatives Properties 


Be f and g derivable in k and ‘c’ a constant. So we 
have to: 


i +g 09-f(9-g (k) 
i. — (f*g)(k-f(k*g(-f(9»eg (k) 
i. — [f (dk) g (À)] 2 LP (k) g (k - f (l9 g' (9] — 
[ (g (k))’], as long as g (k) #0 
iv. (cf) (kK)=c*f’ (k 


All demonstrations here omitted can be seen in 
GUIDORIZZI (2001). 


2.3 Integral 


Given a function f defined in a range [a,b] and k is a 


given number. The sum bellow is the Reimam sum: 
n 


> f(c)AX, onde c € [XX] 
i=l (7) 


When doing (máxAx1) — 0, naturally the sum will tend 
to K, then we will have the following: 


Jim Sh, f(e)ax; =K (8) 


max AX;0 


By setting the limit, given e >0 there is 6 > O that it 
only depends on e: 


KER, f(e)AX) -Kl < z (9) 
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Such number K is indicated by: 
pe l 
f, FGddx (10) 


So we refer to the equation below as the defined 
integral from f to [a,b]. 


K= Jim Shi f(c)AX, = fe fG@ddx (1D 


mar AX;i- 


2.3.1] | Integral properties 
Be f, g integral in [a,b], and A a constant: 


i LIGO)  gGOldx = f? f) dx + f° g dx 
i. f AsfGOdx — Af, fGOdx 
ii, If fc) =0inlabl.therefore f. fG)dx = 0 


iv. If C € [a,b] and f integrable in [a,c] and in [c,d], 
therefore: 


f rodas = [mas — f ros 


2.3.2 Fundamental Theorem of Calculus 


Be f integral in [a,b], if F is a primitive of f in [a,b] 
then: 


.b 
| flx)dx = Fla) — F(b) 


a 


(12) 
Such theorem is only used in the case where the 
integral is said to be defined. 
2.3.3 Integração por partes 


Be f and g defined and derivable in a range I. Note 
that: 


(f « g)'(x) =f (x) * g(x) +£(x) a) (13) 
fGo*g G)-(*g'G)-fG)*gG) (14) 
fe FG) = g (x)dx = f, If * g) GO) — FG) « g(x)] dx. (15) 


Jpg (x) fdr = f GO + g(x) — [2 f£ GO * gGO dx (16) 


Doing u=f(x) and v=g(x): 


du = f' (x) dx and dv 2 g' (x) dx 


Therefore, integration by parts is: 


nb b 
| udu = ur -f vdu 
a “E (17) 
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2.4 Numeric Sequence 


A sequence is a function whose domain is the set of 
natural numbers, that is, it is an endless succession of 
numbers. 


For all n € N, traditionally, a sequence is written as: 
Xn = (Xo, Xi, X2, ...,Xi, ...) 


A sequence can converge or diverge. If the sequence 
converges to a real number ‘’k’’ we say that it is 
represented by a limit, that is: 


lim X, =k 


=Á a 
2.4.1] | Propriedade de Sequéncias 


The same operating rules used with function limits are 
valid for number sequences. See section 2.3.1 of this 
chapter. 


2.5 Séries Numéricas 


Given a sequence (Xn) n€N, which may converge or 
not, we say that an infinite series is the sum: 


Xo+ XIX +... Xp .... 


With the compressed writing: 


ne 
n=1 


If the series converges, we say that the limit of its sum 


(18) 


exists and will be represented by a real S number, or 
rather: 


xc (19) 
2.6 Periodic Functions 


We will begin the study of periodic functions, 
something fundamental to continue the understanding 
about the series under study. We say that a function f: IR 
— R is periodic period t if f (x+t) = f (x) for all real x 
(Ávila, 1999). 


The function cos (x) it is periodic of period 2IT 
because: 


Cos (x + 2TT) = cos (x) * cos (2TT) — sen (x) * sen (2TT) 
With cos (217) = 1 and sen (217) = 0, it follows: 
Cos (x + 217) = cos (x) * 1 — sen (x) * 02 cos (x) 


Be f and g functions of given period t, then f + g and 
f * g will also be periodic. 


Demonstration: 


(f + g(x +t) = f (x + t + g(x + t), because any 
functions enjoy the definition of sum of functions. By 
applying the periodic function definition, we have: 
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(f+ g)(x+t =f (x +t) + g(x +t) =f (X) + g(x) = (T+ gx) 


O que nos mostra que a soma de duas funções 
periódicas também e periódica. A outra demonstração será 
feita de modo análogo veja: 


(fg) = (x + 02 f (x + 0g(x + t) = faga) = (fg)(x) 

Which shows us that the product of two periodic 
functions is also periodic. 

If function f is periodic period t, then kf(x) is periodic 
of the same period, where k is a natural other than zero. 

Demonstration: 

Be f [k(x +O] 2 f [(x - O - (x * O * (x -O -...- (x^ 
0] 2»f(x-0-c-f(x-O0-f(xc-t--..f(x-0-2f(x)-«f 
(X) + f(x) +... + f (x) = kf (x) 

Which completes our demonstration. 


If fis differentiable and periodic from period t, we will 
show that the derived function f is also periodic for the 
same period. 


Demonstration: 
Note that: 
fies a Cao Ra fa (Co PROAC Ra (Co PP 


Be f: R > R a periodic function of period t, and 
integral in any range, we will show that: 


f - E f A 


Where a it is any real number fixed. 


(20) 


Demonstration 


By the fundamental theorem of calculus, we have the 
following equation: 


nü-t 
e =| f=F+)-F@=F@)-F@ =0 


AS 
gl) = Fla +t) — F(a) = F(a) — Fla) 


Therefore, q (x) = Ü we conclude that @ (x) = Cit 
is constant and q? (a) = qu). 


So: 


f a E f : 


2.7 Even and Odd functions 


(21) 


A function f: R > R is even when f (x) = f (—x), for all 
real and odd x when f (x) = -f (—x), for all real x. 
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For example, the function f (x) = x? is even, because f 
(=x) = Cx? = x? = f(x). 


The function f (x) = x? is odd, because f (-x) = (-x)? = - 
x 3= f(x). 


2.7.1 | Propositions 


Now let's study some properties of even and odd 
functions. They are: 


1. The sum of two even functions is an even 
function, the sum of two odd functions 1s odd; 


il. The product of two real even functions is an even 
function, and the product of two odd functions is 
an even; 


iii. The product of a real function odd by an even real 
function is odd; 


iv. Be f a real even function. The defined integral of f 
in [-1,1] is twice the defined integral of f in [0,1]; 


V. Be f an odd and integrable real function, then the 
defined integral of f in [-1,1] is null; 


vi. Every real function f (x) can be decomposed into 
the sum: 


f (x) = fp (x) + fi (x) 


Where, fp (x) is an even function and fi (x) is an odd 
function, defined respectively by: 


fp (x) = (f€ (x) + f (x) )/2 

fi (x) = (f(x) - f (-x) )/ 2 
Demonstration: 
Let's now demonstrate each item mentioned earlier. 
i. Be q:R—R and o:R—R even functions. 


(P+O)(-X) = QU) + o(-x) = q(x) + o(x) = (9 + o x) 


il. Be q:R—R and o:RR even functions. 


(P*O)(X) = P(X) + o(x) = q(-x)*o(-x) = (9*o)(-x) 


ili. Be q: RR an even function and œ:R—>R an odd 
function. 


(P*O)(X) = P(X) + ox) = q(-x)*o(-x) = (o*o)(-x) 
iv. As we know we can decomminate the integration 
limits of any integral defined in a finite range, as 
follows: 


| f x | f (x) dx = [7 Go) dx + f fl) dx 


Doing x = -y, you have: 
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D i D i 
2. | P | III. FOURIER SERIES 
ES fl—ydy + | f(x)dx = | fldx + | f(x) dx 
= D -l D A Fourier series follows the form: 
= +a, cos(x) + b,sen(x) + a. cos(2x) + b, sen(2x) + -- 


Because f(x) is even. Like this: 


Or otherwise: 


D I I I i 

fO)dy+ | fG)dx- | f@ddx+] f@dx=2] fla: - 

[rers [room [rn E Ra cien 
l nzl 


(22) 
Therefore: It's called the Fourier series, with n=1,2,3,4,..., where: 
al el 1 r= | | | 
| f(x)dx- 2 | fGdx a, = -f f(x) cos(nx) dx 
-i p Tlm (23) 
V. The demonstration will be made in a manner i;* ._. — 
analogous to item (iv): b, =z] fO) senl(nx)dx 
i | ms (24) 


: They are called Fourier coefficients, of the 
| f= | Hd | Hds + | flxddx trigonometric series which will be explained in more 
_I -I -i D detail. As we know a series can diverge or converge. If the 
series presented earlier converge sum, it will be a periodic 

function f(x) of period 27, because sen (nx) and cos (nx), 


Doing x = -y, you have: are periodic functions of period 2z. So that f (x+t) = f (x). 


3.] Determination of the trigonometric coefficients of 


-[ fl—yay+ | fixddx= | flyldy «[ Flxde= — [ f (dx + | FG dx roana 
l ü l ü " Be FG) = = 4 eae A cosínx ) + b sen(nx)]. 


According to BUTKOV (1988) we have: 





Therefore: 
I 
| r@dax=o U 
-i üg = — | flx)dx 
Tim (25) 
vi. Note that 1f f(x) is pair, you have: 
iy* . 
"——— r — C a, = -f f(x)cos(nx )dx 
Fi) = jie eee torto Tim (26) 
= = - "P" s 1 m | 
Peca e: fp (x) = f (x) = f (-x), by definition. If f(x) is E e | f(x) sen(nx)dx 
odd, it has — if: Te (27) 
3.1.1 Proposition 
Fa 2f(x) | fo» + fl) At) — f(-x) Tl) — fi-x) 
fix) = 2 |. 2 B 2 i 2 To demonstrate the Fourier coefficients we will make 


use of the following proposition. Suppose the functions are 


integrable in a range I where the series is evenly fit. So: 
Because, fi (x) = f(x)= -f (-x), by definition. From this, 


we conclude that: 


fad + f(a) = BEPC) , Ala + fs) _ f) fts) fled fos) f; X a 2 Alas ve 
fp) + fé +f) —-f(-x) f@)+f) 2ff) 0, 
ee ee Demonstration: 

Note that: 


Therefore: 


fd = pl) + flo) 
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aor . 
fi » «e dx= f,[f GO + f. GO ++ f GOldx | f@ax= AT €» a; =) reas = 


b È fn e dx = ff, GOodx f, £ GO dx +--+ f f dx 
n=1 


To determine Fourier's coefficients a, and b,, we will 
use the following auxiliary integrals, i.e., orthogonally 


3 e 3 2 relationships. If n and k are whole and if nZk, you have: 
E > ha e dx= > [f fold 
n=l n=1 
if 
This demonstrates the theorem mentioned above. Now, | cos(nx) cos(nk )dx = 0 
let's deduce Fourier's coefficients. When a function is -r 


written in the form of a Fourier series, we have the " 
following: | cos(nx) sen(nk )dx = 0 
-E 


00009 X | 
fo = Ej t >. La, cos(nx) + b,sen(nx)] x 
a (29) | sen(nx) sen(nk jdx = 0 
= IE 
The Fourier series can be integrated term-to-term in the 
range (-7.7) due to preposition 3.1.1. So we have: Demonstration: 


We will show the results presented by the auxiliary 
integrals. 


[. fodx- [. E | +f Y tog cost + b,sen(nx)] dx 


(30) E | E 1 i : 
cos(nx) cos(nk)dx = 3 [cos(n +k) + cos(n — k)]dx 
= -E 
Now, let's calculate each integral separately. 1 r7 | 1 f7 
= | cos (n + k) dx + — | cos(n — J dx 
2J x SIS ` 
E a, E a, 
— dx = = | dx==a,m | du 
-r é é dm 2 u = (n +k)x >- = (n +k) > du— (n+ k)dx 
dx 
| | aQ,cos(nx)dx = = | cos(ujdu = — [sen(nm) — sen(—n7)] and 
| | aw | | ! 
w= n- k)x > Ta = (n — k) 2 dw = (n — k)dx 
E — By em 
| a, cos(nx) dx = = [senlnm) — sen(—nz)] = 0 Therefore: 
= 
1 [. 1 it 
z E | —— | cos(u) du + ———— xl cos(w) dw 
| b,sen(nx)dx = E | sen(u)du = E [cos(nz) — cos(—nz)] 2(n t k) 2(n — k) w) 
e ndo | 


E | b, | | Integrating, we have: 
| a, cos(nx) dx =—l[cos(nm) — cos(nm)] = 0 
= n 


Going back to: sen(w)] ^. = 


1 T 1 
2(n+ k) sentuM y 2(n— k) 


| | fGddx= | . dx + | . >. La, cos(nx) + b, sen(nx)ldx 1 T 1 T 
n=1 anri + Km E ao Kil = 


We have: [sen(n + Km] — [sen(n + kiiri] — 


1 1 
cine E) 2(n4 E) 
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1 
= py —— benin — kim] — TEST i [sen(n — k)(—7)]] = > x 20-4 kj [7 cos (u)]*. fuo [- cos(w)] 7, = 
1 1 = 1 ee 
2 K) ziii o 2(n+k) dba CA — 2(n4 Kk) [~ cosin + kje] a ta p (n — kK)x], = 
1 : i | 
f0+o0}=0 1 1 
x - 2(n— k) Rn EHI E 2(n —k) cos[(n — kY—m)]] = 
it = mr) [cos (n + kJ] — 20s i sl t k)y(m]- 
sj cos[(m + kjx] dx -5| sen[(n — k)x] dx = 
P als EU mt 
d 2(n— k) 2(n —k) 
u=(n+k)x > = = (n+ k) > du(n+ k)dx — 
w —(n—k)r2 A = (n - k) 2 dw(n — k)dx Now if nz*k, we have: 
IE | i sen(u) du + na | pm dw | Ze nx )dx 
E k) pues (wl, o = E) [-cos(w)] ^. ES | cos(nx) sen (nx )dx 
1 T 
2n k) [-cos(n + kJ ^ gs T: -E z = [-cos(n — 12 hd | pee ey ae 
n (32) 
1 1 
G (n+ k) Ci m 24 p) ela + k)C79)] — Deducing equation (32) we have: 
= E) [cos(n — kJ] — TRE : 2m-] mh- k)(—0]) = | cos? (nx)dx = [5 = [1 + cos(2nx)]dx 
= k 1 " Ty af" | 
iz (n4 Kj [cosin + k)r] — 24k) [cos (n + k)(2)] — | 5 E | costans) ds 
l 
ak [cos(n — kJm] — uL o [cos(n — Km) 5x +_sen(2ns)]", = E (sx = sen(2na)  — sen(2n(m))) 


1 1 1 
- - - 5 (m +7) + Iz sen(2nz) + 24 79" (2n(n))| 
| sen(nx)sen(nk)dx — | sen[(n + k)x] dx + | sen[(n — k)x]dx 1 1 
ad : lá = 5 (20) + [= sen(2nz) + = sen(2n) | 
u-(n-k)r5 a (n+ k) > du =(n+k)dx , 
T+? msn (2nm ) 
w —í(n-—k)x =" — = (n — k) 2 dw = (n — k)dx 

1 

m +— sen(2nm) 
L| m 2n 
| sen(u)du + | sen(w)dw = 

üi a With O = sem (27) = sen (2nz) and z + (1/2n) * 0 = m, 
therefore: 


1 * | 1 m | 7 
2(n+k) +k) J sence du + 2(n—k) | sem (ur) dw = 
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| cos^(nx)dx = m 
= (33) 


Deducing the equation (33) we have: 
E 
| cos(nx) sen(nx) dx = 
= 


d 
u = senny] 2 E = —ncos(nx)2 du = —ncos(nxjdx 


A. coli = = — [-sen(2nz) — sen(2n(m) | sen” (made 5 J udu = 
ue. sen^(nx) E sen i} sen“(n(—n)) 
-znr = = an M | 2n 2n 

sen lnx) sen (nry 
i —— +t EO 
2n 2n 


Therefore: 


| | cosinx) sen(nx) dx = 0 
C» (34) 


Deducing the equation (34) we have: 


| sen^(nx) dx = | 5 —[1 — cos(2nx)]dx 


*1 * cos(2nx) 
| = dx + | e dx 
| 2 TEE 


E o x 1 =e, Da fr 
: {in — (—-2)] — + E [-sen(2nm) + sen(2n (—1))] 


1 
—-?r-+— 


: — [-sen(2nz) — sen(2n(m)) 


1 , | 
T + in [—2sen(2nm) 


1 
T — 2n 2sen(2n(x) 


As O=sen(21) = 
therefore: 


sem (nx) e m - (1/2n) * O= m, 
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a 
| sen-(nx)dx =a 
-E 


Now let's determine the coefficients an and bn. Suppose 
f(x) is periodic, can be represented in the form of Fourier. 


fix) = = + 3: La, cos(nx) + b, sen(nx)] 
n=1 


Multiplying both members by cos (kx), with k 0, we 
have: 


Flx)cos(kx) = = cos(ks) + o cos(nx)cos(kx) + b,sen( nx)cos(nk)] 


diam 


By Proposition (3.1.1), we have that the Fourier 
trigonometric series can be integrated term to term. 


[. flx)cos(kx) dx = j = cosi) dx + 


"Y La, cos(nx) cos(kx) + b sen (nx) cos(nk)]dx 


E y=1 
With n = k: 


E Flx)cos(kx)dx = a = cos kx)dx + ` T La, cos iny] + by sent nx)cos(nk)] dx 


| fx)eos(kx) dx = | | cos^(nx) dx = am 


-E 
Therefore: 


üp = - | f(x) cos(kx) dx 
diis (35) 


Now, multiplying the two members of the equality by 
without (nx) and integrating again over the range (-z, 7), 
one obtains: 


flx)senlkx)dx = > sen(kx) dx + 
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2. La, cos(nx )sen (kx) + b, sen(nx)sen(nk )]dx 


nzl 


| f(c)sen(kx) dx = b; | | sen^(kx) dx = bpm 


1 r" 
ES | f(x) sen(kx ) dx 
Tlm 
3.2 The Fourier theorem 


If an f(x) function is periodic or not written in the form 
of Fourier then it converges evenly to the midpoint of the 
function at each point. That is: 


fx) = > + >. La, cos(nx) + b, sen(nx)] 
n=1 (3 6) 
Where: 


1 E 
Ep = =| reos 


(37) 

a, = | f (x)cos(nx dx 
TOS (38) 

= z | | f (x)sen( nx )dx 
T-n (39) 


3.2.1 Examples of functions developed in the 
form of Fourier 


Any function can be expressed in the form of Fourier. 
If the given function is periodic we have nothing to say but 
if it is not we will have to define it in a certain subset of 
your domain. Here are some examples of Fourier series 
development. 


Example 1: A periodic 27 period f(x) function is set as 
follows: 


fF) = xnxx 


In this way, we have: 
a 2 = ES =- — — | — — 0 g = ü 
üg = [. xdx = E E. zi E 
Therefore, ao=0. Calculating ak, we obtain: 
1 f" 
Gp =— | xcos (kx )dx 
T — H 
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By doing the part-by-part integration technique: 


ur — | rau 


au 
u =y > — = 15> du = dx 
dx 
dv = cos(kx) dx =x => | dr = | cos(kx) dx 


1 
v= Que (kx) 


-—. lir m | q senlkx) dx = 


xsenlk x) ar 


-——. Vx - [ eos], = 0 = 


Calculating bx: 


"EHE cm = [ coslix)a 
Rl ii = Est) cos. x= 


2 
r_ypyeti— 
(neo 


na alsen sen(2x) | sen(3x) sen(4x) | 
fe =2| RR = =  « 














ad = ke 1 Sen(kx) 
pe 


(kx) 
feo = c ps—— 
(41) 


Example 2: Be f(x) = cos(3t). Calculate the coefficients 
do, dn, bn. 


Calculating ao: 


1 
Ep = =| costa! >a,= xj. cos(w)dw "m 


Gp = — |  cosQw)dw]*, > asy = — [-sen(3r) — sen(—3n)] 


1 
a, = qm senlar) +sen(3m)]=0 


Because the sin(x) function is odd i.e.; -sin(x) = sin(-x). 


Calculating an: 
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1 m 
a, — — | cos(3t) cos(nt)dt 
ae (43) 


a, = -f E [cos(3 — n) t + cos(3 + n)]tdt 


1 ;* 1 [* 
i-o | . cos(3 — n) tdt + 2a | m (3 +n) tdt 


sen(3— n)t]*. — sen(3 + nt] y 


-1 
"n = 3n(3 — n) 2n(3 +n) 


Applying the integration limits, you have: 


— —— [senl G — n)z] — sen[(3 — (=| = 


-E53 = 
— — ZG ru [senIG + mz] — senl@ + n)(-m]] 
—1 
E = — |sen[(3 —n) z] + sen[(3 — nrl] — 
—1 
zz G ru [senlG + me] + sen[G + n)an] 
a mE 1 Ss 
a, = a — sen[(3 — ma -p sen [3 + nal 
Therefore: 
ue nz(3 — n)(3+n) 
Calculating bn: 


1;* 
ba =— | cos(3t) sen(nt)dt 
m^ (44) 


1 f" | | | 
b, = 537 | [sen(3 —njt + sen(3 +n) t]dt 
2m Jo. 
1 r” 1 [E | 
b, = — | sen(3 — n) tdt + — | sen(3 + n)tdt 
2m -E 2m -F 


1 1 
b, = ma- my (3— nit t]*. + Zn +) ^ (3+) Jr. T 


1 1 
ba = mGiy 518 + nlr] -— In(3 +n) 
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Jum B-n) cos[(3 — n) x] — re (3 — n)l—m) 


Because the function cos(x) is even, that 1s, cos (x) = 
cos (-x). We conclude that b, = 0. 


Therefore the Fourier series representing the given 
function will be: 


2nsen (11) 


fl) = neuem cos(nt) 
i n (3 — n)(3 t n) (45) 


IV. FOURIER SERIES APPLICATION 


In order to present the applications of the Fourier 
series, we first have to analyse graphically the functions 
even and odd and the ones that are neither even nor odd. 


4.] Graphical Examples 


To start, we will analyze three examples: the first 
includes an even function, the second an odd function, and 
the third a function neither even nor odd. 


Example 1: Be the par and periodic function of period 
2T: 


{rte =ft,-n<tin 
fl + 2m) = Fl) 


We demonstrate graphically that this function can be 
synthesized in a sum of cosines of the type: 


m i 1 1 1 1 
Fiti=— — 4 (cost — z cos2t —cosst — —cos4t+— cosbt— ^ 
f (E) 3 4 9 16 25 


Each cosine term in the series is called harmonic. The 
first harmonic is called fundamental harmonic. Therefore, 
cos(t) is the fundamental harmonic; cos(2t), the second 
harmonic; cos(3t), the third harmonic and so on. It should 
be noted that the amplitude of each harmonic is equal to 
(—1)n * 4/n2. So the amplitude of the twentieth harmonic 
will be 1/100. Also note that the angular frequency, 
defined by the number of period per unit of time and 
expressed in radian per second, of each harmonic grows 
linearly with the order of the harmonics. 


Figure | illustrates the approximation of a par function 
by the sum of cosines. 
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Fig. 1 approximation of a par function by the sum of 
cosines 


In Figure 1, graph (a) represents the function f(t); graph 

(b) shows the first two terms of the sum overlapped with 

the graph of function f (t): 

— — 4cos (t) 

3 (46) 

In (c), the graphs of the function of the second 

harmonic, cos (2t); in (d) the graph of the first three terms 
overlapped with the graph of the f(t) function: 


Es 


— — 4 cos(t) + cos(2t) 
: (47) 
The graph (e) plots the function and the third harmonic, 
and finally in (f) the graph of the first four terms overlaps 
with that of the function. 


Examining these graphs, it is observed that with only 
three harmonics one already has a reasonable 
approximation of function. Note that the amplitude of each 
additional harmonic decreases while the frequency 
gradually grows. Thus, the small details in the function 
graphics are filled in by the higher-order harmonics. 


Intuitively, we can say that as the number of harmonics 
grows, the approximation of the function, by the sum, 
becomes more and more perfect. 


As it is well remembered by SILVA and 
PINEHIRO (2012), the analysis of harmonic signals using 
the time domain in most cases is of complex resolution and 
with many integrals. An easier way to analyse signals is by 
using the frequency domain using the Fourier Transform. 
The Fourier Transform (equation 48) is a useful tool, 
because through it differential and integral equations are 
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reduced to simple algebraic equations. Its use is of great 
importance to know analysis of the total energy of the time 
series. 


fla) = | f (De at 

tem (48) 
Applying the Fourier Transform, we have a graph usually 
called the Energy Spectrum or Frequency with on the axis 
of the abscesses having the frequency and axis of the 
ordered energy. Each bar of a graphic is a phasor with 
certain frequency and amplitude. Using the frequency 
spectrum graph, we can represent a non-periodic signal as 
the sum of its phase components (SILVA and 
PINHEIRO, 2012). 


Example 2: The function we just analysed is a periodic 
even function. Now, we're going to investigate an odd 
periodic function. Therefore, the following periodic 
function of period 27: 





| f) =-l,-w<tan 
fO=1 Octeax 
ft + 22) = FO) 


Similarly to the previous example, we draw the 
graphics of the function and the first harmonics. As this 1s 
an odd function let's now use sine harmonics. Such an 
approximation and provided by: 

zx P 1 . 1 l1 ` 
f~ z (sint t „5indt + 5 sinat + mt qose ) 

Note that the amplitudes of the harmonics are 
expressed by 4/(2n + Dx, withn>1. 


The graphic analysis in Figure 2 follows exactly the 
same script used in Figure 1. It is observed that with only 
three harmonics it is not possible to adjust this function as 
well as the previous even function. This means that 
convergence is slower in this case. In other words, much 
more harmonics will be needed to obtain the same degree 
of approximation obtained in the previous case. The reason 
for this is that this function is more uncompromising than 
the previous one because it is discontinued, whereas the 
previous one is a continuous function. 
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Fig. 2: Approximation of an odd function by sum of sines. 


It is still surprising that even though it is a 
discontinuous function, it can be approximated by 
continuous functions (sine) with the degree of precision 
one wants. And only increase the number of harmonics. 


At first mathematicians were very reluctant to accept 
this fact proposed by Fourier. At that time even the 
concept of function was not well defined yet. In fact, much 
of the advance of mathematical analysis in the 19th and 
early 20th centuries is due to this problem of convergence 
in the Fourier series. 


Example 3: In the last two examples we have 
graphically seen that a pair function can be approximated 
by cosines and an odd function by sines. What if the 
function is neither even nor odd? We will soon see that 
any function (periodic or not) can be decomposed into an 
even component and an odd component. Therefore, a 
periodic function neither even nor odd can be 
approximated by cosines (even component adjustment) 
and sines (odd component adjustment) simultaneously. Be 
therefore the periodic function of period 27: 


f(t + 2m) = f(e) 


Similarly to the two previous examples, we show 
graphically that this function is approximated by a sum of 
sines and cosines of the type: 


1 1 l 
3 = TE (cost — msint) + E (cos2t — 2msin2t) 


f(t)-2sinh1l 


3 (cos4t — 4nsin4t)—... 


1 . 1 
— FERE (cos3t — 3msin3t) + LE 


1 + 
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With harmonic amplitude (cos(nt)-nz.sin(nt)) equal to 
2sinh(1)/(1+n?.77). 


A 
`~ 
| or a 
vs 


" p" -= 
= 
Du 


Tur 
29 
>. 
o 
mns E 
-n-ae 
Dyji 


+ 
—— - 


, i ! | Py + 4 
pue? i i ' 
| fi 
| "M 
e "E d ie” ef Je M mr 


UT m dE. di LM TL a N 
- - LE 


Fig. 3 Approximation of a function (neither even nor odd) 
by sum o cosines and sines. 


As in the two previous examples, it is observed in 
Figure 3 that as the number of harmonics increases, the 
adjustment between the function graph and the sum of 
harmonics gradually improves. 


A very important detail that should be highlighted in 
the three examples we have just seen is that in the first and 
third examples a constant term appears in front of the first 
harmonic: in the first example, the value 72/3 is obtained, 
and in the third example sinh(1). In the second example, 
apparently this constant term does not exist. The point of 
this is that if the function is anti-symmetrical with respect 
to the axis of the abscissas, the constant term is identically 
zero. This is the case with the function of the second 
example. The functions of the first and third examples do 
not have this type of symmetry, and therefore the constant 
term is different from zero. 


The constant term shifts the function vertically up or 
down. It is called the DC component of the signal, what 
stands for direct current. When you add continue current to 
an alternating current, you get a vertical offset of the 
alternating current. Hence the reason for the abbreviation 
DC. We'll get back to that later. 


In the last example, the fact that the function was 
unfolded in an even component and another odd 
component was fundamental to be able to simultaneously 
use the sine and cosine in the approximation of the 
function. This is a general rule, that is, any function can 
effectively be decomposed into two components, one even 
and one odd. 
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4.2 Application in Electric Current 


An alternating current i(t) = A.sen(x) has passed through a 
complete wave rectifier, which transmits the absolute 
(instantaneous) value of the current. We will show that the 
Fourier series that express is given by HALLIDAY (2003): 


cos(nx ) 





á n*—1 
NIA. (49) 


it it 


If it were a half-wave rectifier it would be: 


cos(nx + x) 
| n^—n 
N= ld... (50) 
However, first we must define the function that will be 
written under function of Fourier: 
ff (+) = A.sen(x), 
f(t + 2m) = f(t), 


QO tion, 
gets 20 


Calculating the coefficient ao: 


2H 


üg = sf Asen(x) dx — -[ ane 


T E 


A o A 2 
ay — zv cos(x)]E + — coslx)]ç" 
ay — E 


Calculating the coefficient bn: 


1 E | | l 2H = | | 
b,-—-— | Asen*(nx) dx — — | Asen(x) sen(nx) dx 
Ty T 


TER 


For n=1, you have: 


1 [* 1p 
b,-—-— | Asen^(x)dx — — | Asen? (x )dx 
| 7T ly 


nad 
b= = I (1 — cos(x))dx — = J | (1 + cos(2x))dx 


A__ sen(ix) _ a 
b= 5,10 t — 5 e —Xx|r — 


sen(Zx) m 
EL. 


A 
b, = —tr-2 m) = 0 
i LU mam) 


If we continue, we will see that all terms in sine are 
zero. Calculating the coefficient b;: 
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1 E | | 1 ZR E 
b, =— | Asen(x) cos(nx) dx — — | Asen(x)cos(nx )dx 
| T 


p DOTT 


Note that it is enough to solve only one of the above 
integrals, because the only difference between them are the 
limits of integration. We will call the first integral (I) and 
resolve through part-by-party integration. 


= - | Asen(x) cosínx) dx = 
L sen (x)sen (nx) — | | Asen nx) cos(x) dx — 


cos(x) cos(nx) I 


—sen (x) sen(nx) Gs 
" T n 
1, 1 — PC 
no. n = 


Applying the integration limits we will have: 











p 14 - ! 1 - ! 
I= 21 © sen(x)sen(nx) — E COS (x) cos(nx)]Z 

n l. 1 

i= — ——(-1)"—- +] 

nº — IL nº n 

n^ 
i= E À ia 1 

n*^—1 pan | 


Where, n=1,2,3,4, ... . So the result of the integral is: 


1 ;* no M —2 
Ias | Asen(x) cos(nx) dx = ————— 
dy mn“ — 1) 


Applying the integration limits of the second integral 
and summing both we have: 
—4A 
WO m(n* — 1) 


Where, n2, 4, 6, 8, ... . Therefore the function 


representing the function will be: 





xA | (51) 
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The function written in the form of Fourier will be: 
rf.) — pdsentxlse üzxem 
f (x ) m | Ü.zenzrzzm 


Calculating the coefficient ao: 


1 FF | 1 FF 
a, = -f Asen (x) dx + -f Od x 
Tdo Tdo 


Ar | 
üg = E | sen(x)dx = 
7 Jy 


—Acos(x) | 2A 


it T 


Calculating the coefficient by: 
1 nit 
ba = — | Asen(x)sen (nx)dx 
7 2 
Where n=1, you have: 


bi = =| ser (x)dx = zl ee lide — - 4 cmm) 


ü m 





ba | e 


Therefore, the first term in sine will be: 


; sen(x) 


If we continue we will see that all terms in sine will 
cancel out, except the first. 


Calculation of an: 
1 ni 
a, — — | Asen (x) cos(nx)dx 
TE. D 


By the previous demonstration we have the result of 
this integral. 


— 2 


Doing n=k+1 in ay: 


—2 


Therefore: 


WWW .ijaers.com 


International Journal of Advanced Engineering Research and Science, 8(4)-2021 


PES n (nx) — =) E = = 
2 2 I E nº +n 
n-1,15.. 

ARAÚJO and MÁRQUEZ (2015) examines a linear 
differential equation that represents a simples electric 
circuit of the type resistor inductor (RL), of which the 
electrical tension is expressed by a function E (t) 
impossible to be resolved analytically. In this particular 
case, the solution have not been expressed by a 
combination of elementary functions due to the nature of 
the function. However, by the use of the Fourier 
Transform and complex integration, it was possible to 
obtain an explicit expression involving the model 
parameters presented earlier in this paper, to represent the 
intensity of the electrical current of the examined circuit as 
a function of time. 


4.3 Application in square oscillations 


Another Fourier application to analyse is the square 
oscillation that can occur in an electronic circuit 
designated to guide the rise of the pulses. Suppose the 
oscillation is defined by: 


TÊ Ü0,se—mrzx-zm 
f (x) m seeing 


Using the equations for the calculation of Fourier 
coefficients, we have: 


1 nit 
a, — — | hcos(nt)dt = O 
Tdo 


h1 
? mn 

If n is odd: 
b, = 0 


If n is even, the result is: 


d h 2h pere E sen(3x) sen Ex) E «| 
2 1 3 3 

Note that all terms in cosine disappear in the range (- 
n.T), that the terms that will represent the interpolated 
function will all be sine, and that this type of application is 
related to the following reasoning: As our electric current 
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is of the sine type, but there are machines that do not work 
with this current (a good example of this are the pulse 
motors), the current arrives at the machine with a sine 
behaviour and the rectifier transforms it into the current 
that the machine needs. This transformation was described 
above through the Fourier series. 


A Fourier series can also be represented as a spectrum 
of frequency. In a graphic, amplitude vs. frequency, as 
proposed in Figure 4 by ANDRADE (2003), the relative 
amplitude of every frequency of square wave are 
represented by simple lines. 


ja) (b) 
Espectro 
de Frequência 





Fig. 4 (a) square wave representation followed by its 
frequency spectrum (b) (ANDRADE, 2003). 


It is also demonstrated by Andrade (2003), how simple 
it is to explore theoretical aspects of physics and 
mathematics utilizing oscillating circuits in laboratory. In 
his experiment, it was obtained the Fourier coefficients 
applying to the circuit a tension with the shape of a square 
wave and frequency wo. Although a square wave 1s applied 
in the circuit, the electric current circulating is sinusoidal 
with maximum intensity in this shape. Therefore, it is 
conclude that a square wave is composed by an infinite 
series of harmonic waves with determinate frequencies and 
in consequence is possible to create a Fourier spectrum for 
the wave. 


For the numeric solution of different equations in the 
context of analysis of oscillators amplitude, whenever is 
applied the Fourier transform the results a simple and 
precise (CARLOS EDUARDO, 2018). In the case of 
resonance the Fourier transform is again an important tool. 


V. CONCLUSION 


We hope that the reader has understood the process of 
expanding a defined function in a range in an even or odd 
periodic function. The choice of the series in terms of sines 
(odd expansion) or cosines (even expansion) depends on 
the nature of the problem that one wants to solve. 


It should be in mind that the periodic expansion of a 
function restricted to a range serves only as a support to 
use the Fourier series technique; and that once the 
expansion is done, the values outside the definition range 
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of the function become irrelevant, and therefore 
disposable. Only the approximation of the function, by the 
Fourier series, within the definition range of the function is 
of interest, the rest, I repeat, is disposable. It is good to 
know that, in practice, non-periodic functions defined in a 
range are more important and frequent than the periodic 
functions themselves. 


So far, the impression one has is that any periodic 
function allows to be approximated by a series of Fourier. 
Is this true? The definitive answer to this question is subtle 
and is outside the objectives of this work. However, we 
can ensure that all periodic functions found here, in 
particular in electromagnetism can be developed by 
Fourier series. Technically, they meet  Direchlet's 
conditions, which are sufficient conditions for a periodic 
function to be expanded in Fourier series. 


Periodicity is one of the prerequisites for Fourier serial 
development. This leads us to ask the following question: 
If the function is non-periodic, is it still possible to speak 
in Fourier series? Roughly speaking, the answer is 
affirmative. There are two possibilities to be analysed: (a) 
functions restricted to a finite range of the line and (b) not- 
periodic functions defined in any straight or semi-straight. 
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